LONG TIME SEMICLASSICAL APPROXIMATION OF QUANTUM FLOWS: 
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Abstract. Let 7i be a holomorphic Hamiltonian of quadratic growth on IR^", b a holo- 
morphic exponentially localized observable, H, B the corresponding operators on L^(]R"') 
generated by Weyl quantization, and U{t) = expiHt/h. It is proved that the L? norm 
of the difference between the Heisenberg observable Bf = U{t)BU{—t) and its semiclassi- 
cal approximation of order N — 1 is majorized by K N^^"^^^^^ {—Mogh)^ for t G [0,T]\[(h)] 
where T^ih) = — ^ . Choosing a suitable N{h) the error is majorized by Ch^°^^^°^^\ 
< t < I log ^1/ log I log ?i| . (Here C are constants independent of N^h). 



1. Introduction and statement of the results 

Denote O := IR^'" with coordinates (x,^). Let n{x,^) G C°°(0; IR), and btix,^) := 
bo (j)^ = b{4>Y{^ii)) be the time evolution of any bounded observable b{x,^) G C°°(0; IR) 
under the the flow (pl^ : O ^ O generated by the Hamiltonian H. Denote H := Op^ {H) 
and B = Op^ {b) the self-adjoint operators in L^(IR"^) representing the (Weyl) quantization 
of the symbols H.b and let Bt := e'^^^^Be~''^^^^ be the Heisenberg observable, i.e. the 
quantum evolution of the observable B under the unitary group generated by H. 

The question of estimating how long the classical and quantum evolutions stay "close" 
one another or, better, how long the evolution of the quantum observables is determined 
by the corresponding classical one up to a prescribed error vanishing with % is one of the 
oldest problems of semiclassical analysis. According to a well known conjecture going back 
to Chirikov and Zaslavski [Ch,Za], this approximation can be valid on a time interval of 
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maximum duration T = T{%) of order —\og%, called the Ehrenfest time, if the error is 
required to vanish faster than any power of %. 

The origin of this conjecture, formally verified in some instances[Za] can be under- 
stood in the correspondence between symbols h{x,^) (classical observables) and operators 
in Hilbert space B (quantum observables) provided by the Weyl quantization procedure: 

{Bu){x) = J^^^ b e'^^^-y^'^^/''u{y)dyd^, u e 5(1R") (1.1) 

In this framework the problem can be formulated as follows: Bt solves the Heisenberg 
equation of motion 

Bt='-[H,Bt] (1.2) 
If Bt admits a symbol, denoted bt{x,^;h), by (1.2) it fulfills the equation 

bt = {n,bt}M (1.3) 

with the initial condition bo{x, ^;h) — b{x, ^). Here {/, g}M{x, ^) is the Moyal bracket of the 
two observables f,g e C°°(IR^'') 

{f.9}Mix,0--=f*9-9*f (1-4) 

where /#^, the symbol of the operator product FG, is expressed by the composition of the 
symbols / and g: 

ifMix.O = to L/2 / e-'^-^'^y''+'^^^^y^f{x+w,p+i)g{x+T,T+i)dpdTdTdw (1.5) 
(27r/i)'^/^ i]R4" 

{fi g}M admits the following formal expansion in powers of % [Fo,Ro, Vo] : 

{f.9]M{x.i) -{f.9] + ^ E (^1)""^' i^tgD^^g) ■ {dlgD-j) (1.6) 

|a+/3|=i>l 

{a = (cKi, . . . , an) is a multi-index, and \a\ := ai + . . . + a^, analogous definitions for 
and Dx := —ihdx). By (1.6) the differential equation (1.3) can be recursively solved in the 
space of the formal power series in h (for details see [Ro], Chapt.IV.lO). The result, known 
as the semiclassical Egorov theorem, is the formal semiclassical expansion of the symbol bf. 

oo 

btix,^;h) ~ ibo<p-){x,0 + Y.^jix,C;t)n^ (1.7) 



Here the term of order zero in h, by definition the principal symbol of B^, is just the evolution 
of the observable b along the Hamiltonian flow generated by H, i.e. the solution of the 
Liouville equation bt = {H, bt}, and 

b^{x,C;t) = -i f (1 - (-l)'"+^')r(«,/3) {d^nD%)oct>'^_^{x,OdT (1.8) 

0<(<j-l 

The higher order terms bj{x,^;t) are thus completely determined by the classical evolution 
but have a polynomial dependence on the derivatives of the flow 4)^{x,^) with respect to 
the inital conditions (x,^) up to order j — 1. If, as it happens in general, there are initial 
conditions (x, ^) generating a flow with positive Lyapunov exponents, the difference between 

N 

the symbol bt{x,^;?i) of Bt and any prescribed approximation (6o(/)^)(x, 0^"' 

is expected to increase exponentially in time: hence it can vanish as h ^ only for a time 
interval not exceeding — log^. Put in a different way: the non-local nature of quantum 
mechanics, embodied in the symbol expansion (1.7), (1.8), can be dominated by its local 
approximation, the principal symbol bo (f)^(x,^), only if the the remainder is small. This 
can be obtained only within the above time span. 

In this paper we work out, in the analytic case, the estimates implying the validity of 
the above "Ehrenfest time" for a class of flows somewhat restricted but in a sense natural 
as discussed below. More precisely, for any flxed cr > set \z\ := supl^^l and : = 
G C'^^ : llm^l < cr}. The Hamiltonian H{x,^) = 'H{z), z :— {x,^) is required to fulfill 
the following properties: 
(Al) There exists u > such that Ti is real-holomorphic on Q^. 

(A2) Let JdTi be the symplectic gradient of Ti. Then there are Ai > 0, A2 > 0, a > such 
that \Jdn{z + iy)\ < Ai + A2\z\ Wz,y G IR^'", \y\ < a. Moreover \Jd'^n{z) \ < a on Q^. 

(A3) Denote {'H){k) the Fourier transform of T-i{z). Then there are p > 0, cr > such that 
7i(/ci + ik2) is holomorphic on Qp \ (0, 0); moreover k^'H{k) is holomorphic on Qp and 

\kif\n{ki + ik2)\ < Ce-'^l'^il for j/csl < p 

Remarks. 

1 Under the above assumptions H = Op^ {Ti.) defined by (1.1) is essentially self-adjoint 
in L^(1R"^). By a standard abuse of notation we denote H also its self-adjoint closure. 

2 Within the analyticity and decay assumptions (Al)-(A3), (A2) is the quadratic growth 

condition ensuring the existence of the Fourier integral operator representing the pro- 
iHt 

pagator exp — — [Cha] and thus the existence of the symbol of B^ [Ro] . 
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3 In the phase variables z= (a;, ^) Assumption (A3) means that there are cr, p > such 
that 



sup 



e^l^l < +00. 



(1.9) 



To state the main result of the paper we need some further notation. For b as above set: 



Anb:=^[{b,n}-{b,n}j^] . 

and define recursively the two sequences r^, : /c > 1 in the following way 



(1.10) 



(1.11) 



b\ 



fdr, j 
Jo Jo 



t — Tl — T2 l-t — Tk-l 

dT2 I dTs... 

Jo Jo 



Moreover, let 6 : C be holomorphic. Set: 



t-Tl,Tl,T2,...,Tk-l ITk 



o (f)'^'';bo :=bo 



\b\^^p := sup \b{x + iy)\e''^ 



(1.12) 



Denote Acr,p the set of all functions / holomorphic on Qa;p such that |/|^ ^ < +00. Then: 

Theorem 1.1. Let there exist a > 0, p > and < B < +00 such that \b\^ p < Then: 
(1) The operators := Op^(6* ) are continuous in and the Heisenberg operator B^ = 
U{t)BU{—t) admits the expansion 

N 

B, = Y,Bp'^+h'^''+'^Sl,, 

3=0 



where 



qt 

On 



:= / dn / dT2 / drs... / dTkU{Tk)Op'" h 

Jo Jo Jo Jo ^ 



t— Ti,ri,T2,...,rfe. 



ui-n) 



(2) There are positive constants E, F independent of j, N and % such that for all j > 1, 
N > 2, t > the following estimates hold: 



I R*ll 



exp a 1 



6n+3 



N(N-l) 
exp I a 1 



6n+3 



(1.13) 



(1.14) 



Remark. The holomorphy assumptions are needed to control the remainder to order 
for all N. If we limit ourselves to N = 1 more general classes of Hamiltonians and of 
observables can be considered. More precisely, let for instance H{x,^) be a polynomial of 
order 2p such that the subgraph := {{x,$,) e JR^'^\H{x,^) < E} is compact for some E, 
and let b{x, ^) e Cq°(S£). Then (proof in the next section) there are F > and A > such 
that: 

\\B' - ^oIL^^L^ = \\B' - Op^(6o 0*)||^,_^, < Thhe''' . (1.15) 

The symbols 6* and hence the operators are completely determined by the classical flow 
0* via (1.11). The quantum evolution will then stay close to the (semi) classical one as long 
as the error Sj^ stays small. The estimate (1.14) yields indeed, through a straightforward 
computation: 

21 h 

Corollary 1.2. Let Tiv(n) := := Bp\ Then, forO<t< TNih): 



Remarks. 

1 If the Lyapunov numbers are zero for any initial datum (x,^), then we can take a = 
in formula (1-14), and by Assertion 2 of Theorem 1.2 one has 



Bf-B 



N{h) 



o{h^) o<t< TNih) 



-1 ^^-6n-l^-l 



where Tiv(n) := e'^N 

2 Estimates valid for a time interval of duration — Cjvlog/i for Hamiltonians admitting 
polynomial growth of any order (but without control of the constant Cn), have been 
obtained by Combescure and Robert [Co-Rol] in a weaker sense, i.e. comparing classical 
and quantum evolutions along coherent states (according to ideas introduced in [He], 
[BZ] and developed in [Ha], [BIZ], [Co-Ro2]). 

3 The symbol expansion generated by Assertion (1) of Theorem 1.1, namely 

h{x,^;n) = f2bp'^ + o{h'''+') 
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differs from (1.7) in all terms with j > 0. This difference makes the present expansion 
a non formal one, so that its remainder can be estimated. 



4 Finally, let T{h) e C([0, 1]; IR+) be an increasing function such that lini — — - 



= 0, 



and let N{h) 



logh 
T(n)\ 



. Then clearly 



B,-B 



Nih) 



= o{h°°),o<t<T{h). 



To put this result into a more quantitative version, define the function sequence {log^'^^(x)} : 
e IN by log[^^(x) := log(x), log^''\x) log(log[''"^^ (x)). 

Corollary 1.3. For any integer k >1 deGne Nk{?i) := [log'^^^d log(/i)|)]. Then there exist 
positive constants C, % such that, for < h < h one has 



Bt-B, 



Nkih) 



< (7;^iogW(|iog(a)|) 



for 



0<t< 



iog{h)\ 



logW(|log(7i)|) 



Acknowledgments. We thank A.Martinez and D.Robert for reading the paper and sev- 
eral useful remarks. We acknowledge the support of CEREMADE that made possible the 
collaboration leading to this work. 



2. Proofs 

Let 6 be a Weyl symbol of class Eg, and H an admissible semiclassical symbol (For 
these notions, see [Ro], Chapter 2; particular examples are all bounded observables b{x,^) e 
C~(]R^") and the Hamiltonians H e C°°(]R^") of polynomial growth at infinity). Denote 
0* the flow generated by JdH, J the unit 2n x 2n symplectic matrix; let H = Op^ (H) 
be essentially self-adjoint in L^{JEC') and denote also U{t) := exp{itH/h), B := Op'^ (b), 
Bt := U{t)BU{-t), and 

Anb:^^[{b,n}-{b,n}^] . (2.1) 
ft 

Our semiclassical expansion is generated by the following simple remark: 
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Lemma 2.1. The following formula holds 



Bt := Op^ {b o 0^) [ dTU{T)Op'" {r{-^) U{-t) , 

Jo 



where rf := A^^ (6 o 0*) 

Proof. Denote Pt '■= bo 0*. Then: 



I [Op- (A)] = op^ = Op- m,n}) 
= op^ ({A, n} - {A, ^}m) + ^ [^^p" (A) > ^^p" m 



n 



[Op^ m,Op^ {n)] + h'Op^ {r\) 



It follows 



^ [Op- (A) - St] = ^ [Op- (A) - Bt, Op^ m + n^Op- {r\) , 
and by the variation of parameters formula 

Op- (A) - St /i^ _ (^..^ _ ^^^^^ _ ^2.2) 

The assertion is now proved performing the change of variable r = t — s in the integral. □ 
Proof of formula (1.15). Since is a polynomial of degree 2p 



2p 



r\ = An{bo 0*) ^^[{bo H}-{bo H}^] = c^/i^ ' ^ 

" ifei=i 

fc = (fcl,...,fc2„) 



where Ck{x,$,) is a polynomial of degree 2p — Now the smooth functions 9k{x,^) := 

Ck{x,$) — T—r — have compact support in IR^"^ and hence define bounded operators in 
upon Weyl quantization. Denote X{x, ^) the Lyapunov number of the trajectory (p^ with any 
initial datum (x,^) G S^. Since 4)*{x,^) is bounded Vt G IR we have (see e.g.[Ce], 3.12) 



S := supj^ A(a;,^) < +oo. Hence there are 7^ > such that sup^^ 



< 7fee'^*. Since 



— —r — is a polynomial of degree A; in the variables — , s = 1, . . . , /c with coefficients 

oz'^ oz^ 

d\'\b 



depending on 



dz^ 

Tfc ip) > such sup ^ 



, s = 1, . . . , for any fixed q E ^ depending only on n there are 



dz^ 



< Ff.e^'^^'^^Wll < q. Hence by the Calderon-Vaillancourt 



theorem there exists q > such that \\Op^ {9k)\\L^^L^ < F/ec''^'^^*. Inserting this estimate 
in (2.2) we get (1.15) with A = 2pq5, T = Max^Ffc . □ 

Recall now the definition of the sequences '^i.'^i.T2v,-rfc-i ^ 2) and {bl.}f.yQ, bo := b: 



t — Tl,Tl,T2,...,Tk A 



n 



^t-Turi,T2,...,Tk-l ^ 



bi 



:= / dri / dT2 / drs... 
Jo Jo Jo Jo 

Lemma 2.2. Let B] = Op^{b*). Then: 



N 



j=0 



where 
Sn 



:= jy^J^ "'cZrs^ ^'d^^-J^ dTfcC/(Tfc)Op- (r*""^'"^'"^' -'"''-^) C/(-Tfc) 



Proof. Just iterate the proof of lemma 2.1 □ 
Let 6 : ^o- — ^ C be an analytic function; recall the definitions 

\bl^^:= sup \b{x + iy)\eP\-\ . (2.3) 

and A(T,p '■= {b holomorphic in Qcr '■ \b\^ p < +00}. We will estimate the sequence rk in the 
above norm. Clearly we have to estimate the norm of bo (p^ and of A^6. We first prove the 
following 

Lemma 2.3. There exists a positive a such that (j)^ extends to a complex analytic function 



Proof. Denote / := JdTi, and consider, on ^0-5 the system of equations 

z = m . 
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(2.4) 



Writing z = x + iy and / = /i + if 2, one has y = f2ix + iy). Since /2 = on the real axis, 
by assumption A2 one has \f2{x + iy)\ < a\y\. It follows that the inequalities 

\y\<Oi\y\ ^ < |yo|e"l*l (2.5) 

hold. So one has (t)*{Q^^-oc\t\) C Qa- 

Fix t. Given z G Q^^-a\i\ we prove that 0* is analytic at z. By the Cauchy- Kowaleskaya 
theorem (see e.g.[Pe]) there exists a neighbourhood U of z and a time i such that, for any 
|t| < t, (j)'^ is analytic on U. Assume that t is the supremiim of such times (so that (p* is 
not analytic in U). Assume by contradiction t < i. By (2.5) the limit lim^^i (j)'^ (z) exists 
on U. Denote w := lim,-_>f 0'^(z). Again by the Cauchy-Kowakeskaya thoerem there exists 
a neighbourhood V of w and a ti > such that 0"^ is analytic on V for |t| < ti. Assume 
that U is so small that for fixed e small enough one has (j)^~^{]A) C V, then one has 

/+^(W) = 02^(/-(W)) , 

which is analytic since it is the composition of two analytic functions, against the assumption 
that t is the last time of analyticity. □ 

Lemma 2.4. Let b e Aa,p, then, for any t, and for a small enough, one has 6 o 0* e 



Proof. By the above lemma h o (j)^ has the required analyticity properties. Denote pt '■= 



pe 



-a\t\ 



, at '.= ere "1*1, (f\ + i(f2 = <l)*{x + iy), then one has 



60 



sup 



b{(p\x + iy)) eP'\ 



< sup 



using the equation of motion and A2, one has 



at 



\Re{(f) *{(pi + i(p2))\ < \<Pi\e 

which implies the assertion. □ 
We will estimate the norm of A-^ using the Fourier transform. For this reason the 
following lemma is useful 
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Lemma 2.5. One has 



2\" 1 



(2.6) 



Proof. Fix ki = kci where ei is the unit vector of the first axis and k a positive number; 
fix also k2 with |/c2| < p — 5. One has 



(27r)" b{ki + ik2) 



J R2n 7 J^2n 



.-/3|a;|g-«;CTg|fc2||a;| 



2n 



which by definition of \b\^ ^ is the thesis in the particular case just considered. The general 



case can be dealt with in a similar way. 



□ 



Lemma 2.6. Let b e Aa,p with a <v small enough. Then there exists a positive constant 
A such that, \/d < a, 5 < p: 



Proof. To obtain the estimate via the Fourier transform we first recall that 

{6, n}i, {k)^l I b{k - s)n{s) sin ^^~'l^' ds 
nj^2n n/z 

where {kp, kg) A {sp, Sg) := kp ■ Sg — kg ■ Sp, whence 

X — 2 f , / . (k-s)As (k-s)As\, 

Since | sinz — z\ < Ci\z^\ for all z e Qf^, one has, for \lmk\ < p — S, 



Anb{k) 



< 



b{k-s) n{s) 



Ci\k - splspc/s 



\ki - sre 



3„ — (tI/ci— s|„— alsl 



ds , 



(2.7) 
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where use has been made of (2.6) and Assumption A2. Now — s| + |s|| > and 
— s| + |s|| > \ki — s\. Hence (2.7) does not exceed 



J R2« 



which gives 

Using again (2.6) to antitransform the assertion is proved. 



□ 



Lemma 2.7. Assume \b\cr,p < B for some positive B, cr, p. Then, for k >1 and < < t, 
one has 



J^ — Tl,Tl,T2,---,Tk- 



((T-fe5)efc,(p-fcd)efe 



Here the sequence Ck is defined by 



ei := e , := ei exp (—at(k — 1)) , k >2 



(2.8) 



and the sequence Tk by 



1 ■ ^2n(^4n+3 g6n+3 



Tk := Ti 



Aei 



fc-l 



exp I a 



k{k-l) 



6n+3 



(2.9) 



Proof. The expressions of ei and Fi are a direct consequence of lemmas 2.4 and 2.6. By 
induction assume that the estimates of the lemma are true for k we prove them for k + 1. 
By lemmas 2.4 and 2.6 we have 



and therefore 



^^t — Tl,Tl,T2,...,Tk 



^=+1 (fe+l)5)efce-'^^fe,(p-(fc+l)d)efce-''-^fc 



(a-(fc+l)5)efee-'""fe,(p-(fc+l)d)efee-''-^fc 



(efce-"*)6"+3d4n+3^2n 
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This yields Cfe+i = e^e and therefore (2.8); moreover 

A 



whence 



A 



k-l 



6n+3 



(^4n+3^2n 

This proves (2.9) upon insertion of (2.8). This proves the lemma. 
Lemma 2.8. For any N >2 one has 



□ 



Op™ (r 



t — Tl,Tl,T2,...,TN-l 

N 



< 



N{N -1) 



6n+3 



(2.10) 



where E, F are positive constants independent of N. 

Proof. We estimate the l.h.s. of (2.10) by the norm of the Fourier transform of vn- By 



lemma2.5 we have, for k e JR^^, 



^t — Tl,Tl,T2,...,TN-l 
' N 



5-L N 



(p - Ndy^'ej^ 



e'x.p[—{a — N5)eN\k\] , 



and therefore 



^t — Tl,Tl,T2,...,TN-l 
' N 



< 



exp [-(cr - N5)eN\k\] dk 



5-L N 



(p - iVd)2"e^ {a - iV5)2^e 



AT 



Choosing 5 = a/2N and d = p/2N, and inserting the expressions of ejv and the assertion 
is proved because, (see e.g.[Ro], Corollary 11.19) 



Op"" (r 



t— ri,ri,T2,...,Tjv- 
AT 



-') 



< 



U — Tl,Tl,T2,...,TN- 



□ 



Lemma 2.9. For all N >2 and t > the following estimate holds 



exp I a 



N{N -1) 



6n+3 
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Proof. One has 



SnW < 



sup 

0<Tfe<t 



Op^ (r^-^-^-^--^--) II 1^ dr, £ dT2 £ dT3... ^ 



t — Tk-l 



dTk ; 



the norm of the Weyl quantization of tn is estimated by the above lemma. To compute 
the integral it is convenient to make a change of variables introducing the new variables 
si, sn defined by 

Sl = t — Ti , Si + S2 = t — T2 ... Si + S2 + ... + Sn = t — Tn , 



This transforms the integral into 



t f-t-Si rt-(8l+82) /■t-(si+S2 + ---SjV-l) 

dsi / ds2 / dss... / dsN 

^0 Jo Jo 



To see this fact it is enough to remark that 

pt—TN—i r 
I dTN = 

Jo Jo 



t—TN-1 rt — TN-1 rt—{si+S2 + ---SN-l) 

dTN = / dSN — dSN ■, 



(2.11) 



and to iterate the argument. Denote In the integral in (2.11). We claim that /jv(i) = /N\. 
To this end remark that one has 



lN+l{t) = lN{t-s)ds. 

Jo 



Hence the assertion is proved by induction and the result is thus obtained. 
Proof of Theorem 1.1. It is enough to apply Lemma 2.9. 



□ 
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